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Abstract 

Multipersistence homology modules were introduced by G.Carlsson 
and A.Zomorodian \X\ which gave, together with G.Singh [4], an algo- 
rithm to compute their Grobner bases. Although their algorithm has 
polynomial complexity when the chain modules are free, i.e. in the 
one-critical case, it might be exponential in general. We give a new 
presentation of multipersistence homology modules, which allows us to 
design an algorithm to compute their Grobner bases always in polyno- 
mial time by avoiding the mapping telescope. 

Mathematics Subject Classification (2010): 55U99, 55N99, 13P10. 
Keywords: Computational Topology, Commutative Algebra, Per- 
sistent Homology and Grobner Bases. 



1 Introduction 

The theory of persistent homology builds a bridge between computational 
algebraic topology and data analysis using homology as an effective tool to 
associate a computable invariant to a point cloud. Other applications of the 
theory range from coverage problems in sensor networks to complex network 
theory. The main idea is to approximate the point cloud embedded in a 
metric space with an increasing sequence of simplicial complexes (filtration) , 
see [11 [8]. By analysing the persistent, i.e. long living, homological features 
in the filtration, the shape of the point cloud can be inferred. 
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supported by by PRIN 2009 "Spazi di Moduli e Teoria di Lie". This research has been 
partly supported by the TOPDRIM project funded by the Future and Emerging Tech- 
nologies program of the European Commission under Contract IST-318121. 



Multipersistent homology is a generalization of this theory in which the 
homology of a filtration of simplicial complexes indexed by N** (multifiltra- 
tion) is analysed. A multifiltration represents a simplicial complex evolving 
along r different directions. 

Multipersistent homology has been introduced by Carlsson and Zomoro- 
dian in the seminal paper [1], where the authors study the classification 
problem for multipersistence homology modules. From then, several papers 
have considered these modules: the same authors and Singh in ^ devel- 
oped an algorithm for computing their Grobner bases [H |9] and Chacholski 
et al. [6] gave a combinatorial multigraded resolution for multipersistence 
homology modules. The classification problem has been studied in [llj . 

A particular type of multifiltration is highlighted in |3] and called one- 
critical. Intuitively, in the one critical case, simplices enter the filtration in 
one single instance rather than in non comparable states. The algorithm 
in [3] works on one-critical multifiltrations optimally and on general ones, 
reducing them to the one-critical case by using the mapping telescope. 

As observed by the authors in [1] the mapping telescope increases the 
size of the multifiltration exponentially in the worst case. 

Therefore, although the algorithm in |4] has polynomial complexity in 
the one-critical case, the use of the mapping telescope introduces a bot- 
tleneck in the general case. In this paper, we give a new presentation of 
multipersistence modules. By using this, we are able to extend a version 
of the algorithm presented in [3] to all multipersistence modules avoiding 
the use of the mapping telescope and therefore eliminating the bottleneck 
of the previous approach. Our algorithm has polinomial complexity for all 
multifiltrations and in the one-critical case it essentially coincides with the 
one in 

We will now sketch the structure of the paper. Section 2 introduces a 
construction that generates a non one-critical multifiltration from a point 
cloud, motivating our algorithm. In section 3 we lay the theoretical basis: 
we show in fact how the modules of cycles and boundaries are isomorphic 
to submodules of a finitely generated free module. Embedding cycles and 
boundaries in the free module we can adapt the algorithm in [1] to the 
general case. The algorithm presented by Carlsson et al. is described in 
section 4 along with some definitions from commutative algebra. Section 5 
presents our algorithm for the computation of Grobner bases in the general 
case. 
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2 Motivating Construction 



Multifiltrations arise from data sets in many ways. For instance, in analyzing 
the structure of natural images the data is filtered according to density, |5]- 
The following generalization of the Rips Vietoris complex determines a non 
one-critical multifiltration from a data set. 

Let 5 C be a point cloud. Consider the ellipses centered in the points 
of S with semi-axis of lenght a and b respectively. Fixed a direction u S , 
we assume that the semi- axis of length a is parallel to v. Fixing the values 
of a and b, we can build a simplicial complex on S by following the same 
procedure used for the Rips- Vietoris complex. Given n points in S, they 
define a (n — l)-face of the complex if the corresponding ellipses intesect two 
by two. According to this construction, the critical coordinates for (n — 1)- 
faces correspond to couples (a, b) such that n ellipses of semi-axis a and b are 
tangent two by two. As a and b vary in [0, oo) we obtain a non one-critical 
multifiltered complex. Given two points in S there can be non comparable 
values of the semi-axes such that the corresponding ellipses are tangent. 



In fact there are infinitely many (see example 2.1). Note that for a = b 
the construction is the Rips- Vietoris complex and that this multifiltration 
shows if some points are aligned in direction v. 

Example 2.1. Consider the origin and a point on the unitary circon- 
ference centered in the origin {cos{t), sin{t)) t G [0,27r). We assume v is 
the X-axis. For t varying in (0,7r/2), the implicit function between a and b 
in the tangency point of the ellipses centered in the two points is a branch 



of hyperbola, see figure 2. 1 . The hyperbola degenerates continuously to the 
vertical and horizontal axes for t = and Tr/2 respectivelly. The example is 
quite general: the function is symmetric with respect to the x and y axes. 



3 Theoretical setting 

Let us consider N*" equipped with the product order, i.e. v = (vi . . . Vr) ^ 
w = {wi . . . Wr) iff Vi < Wi for all i. We denote hy R = k[xi . . . Xr] the 
polynomial ring in r indeterminates over a field k and x" is the monomial 
x^ •••x^''. Unless otherwise stated, we will call module a i?-module and 
vector space a /c-vector space. 

A simplicial complex X is called multifiltered if we are given a family of 
subcomplexes {X^} with v ^ W ^ such that X^ C X^] whenever v ^ w. 
The family {X^} is called a multifiltration of X. If X is a finite simplicial 
complex, every multifiltration of X stabilizes, i.e. there exists a multi-index 
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Figure 1: Implicit function between the semi-axes a and b in the tangency 
point of the elhpses. 

v' = {v'l . . . v'^) such that Xw-^^^. ~ for ah w = {wi . . .v'^. . . Wr) S N'' and 
1 < i < r, where e^'s are, as usual, the vectors in having one in the i—th. 
entry and zero elsewhere. 

We fix a filtered d-dimensional finite simplicial complex X with multifil- 
tration The vector space Cn{Xy) of n-chains in X^ will be denoted 

by Cn{v). This is the vector space with linear basis B„(f ), the set of n- faces 
in Xv The n- faces of X-^i ~ X will be ordered and denoted by a"^^. . . 
(or simply by cjj when n is fixed) where d„ is the rank of Cn{v')- 

The multipersistence n- chain module (first defined in |4j and described 
in [6]) is the multigraded vector space 

Cn = ®Cn{v) (3.1) 

V 

equipped with the left i?-action defined by the linear maps : Cn(») — >■ 
Cn{» + u) induced by the inclusions X, — )■ X,^u- The set of bases {B„(7;)} 
is closed under the module action. Indeed, x'"B„(») C B„(» + w) for all 

w G 

Definition 3.1. The module Dn := 0^, C„(v' + w). 

Remark 1. The module Dn inherits the grading from Cn in the obvious 
way and it is a free module, isomorphic to R'^"- , with R-basis B„(f;'). All the 
generators have multidegree zero. 



The multipersistence homology modules are the homology modules of the 
multigraded module chain complex (see [W]) 

C: ^ Cd % Cd-i ^ ■ ■ ■ ^ Ci % Co ^ 

obtained as the direct sum of the simplicial chain complexes of for all 
V £ W . As usual, we denote by Zn ^ C„ the kernel of 9„ and with Bn-i ^ 
Cn-i the image of 5„. Since multiplication by monomials is injective in C„, 
the chain complex C, is isomorphic to the shifted chain complex 

x"' C7. : ^ x"' Cd % x""' Cd-i ^ > x""' Ci \ x^' ^ 0. 

Note that in general x'C„ is a proper submodule of D„. It is indeed enough 
to observe that x"Cn{{Q, • • • ,0)) C Cn{v'). 

Since dn{Dn) C Dn-i, as it easy to check, the following is a chain com- 
plex of multigraded modules 

D,:O^Dd% Dd-i ^ yDi%Do^{). 

where the boundaries are obtained by restriction. We see x""' C, as a sub- 
complex of D,. 

Definition 3.2. Let a G B„(i;') be a basis element of Dn, a critical co- 
ordinate for a is a minimal v G N*^ such that there exists t G B„(t') with 
The element r is called a fundamental element associated to 

a. 

In general, critical coordinates and fundamental elements for a G B„(f') 
are not unique. We denote by the set of fundamental elements associated 
to 0" G B„(f') and by the set of corresponding critical coordinates i.e. 
the multi-degrees of elements in J>- The cardinality of these sets is denoted 
by kfj := = \Ca\- The set Tn = Uo-eB„(?;') -^o- generates the module Cn- 

A presentation of C„ is given by the short exact sequence 

where i?*^" is the free module with basis the fundamental elements. We define 
maps dn '■ R'^" — T- Dn-i such that the following diagrams are commutative. 
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(3.2) 



The maps dn are linear maps between free modules and they can thus be 
represented by matrices with coefficients in R. The columns of such matrices 
are the boundaries of n-dimensional fundamental elements, shifted in Dn-i- 
Each column is given by monomial entries and two elements can differ only 
by scalar coefficients, being the columns omogeneous elements of Dn-i- 

Theorem 3.0.1. The modules Zn and Bn can be calculated from matrices: 

Bn ^ x" Bn = Imc?„+i and Zn ^ x" Zn = x" o (^„(ker9„). 

Proof. The modules Zn and Bn are isomorphic to x"' Zn and x"' Bn respec- 
tively, by the injectivity of multiplication by x" . Considering also that di- 
agram |3]2] is commutative and that <f>n is surjective, the following equalities 
are satisfied. 

x"'Bn = X"' o dn{Cn) = 5„ o gf' {Cn) = MR''"). 

x" Zn = x" ker(9„) = x" ker(x^ o dn) = x" o 99„(ker5n) 

□ 

Example 3.1. Consider the multifiltered finite simplicial complex with four 
vertices, five edges and one 2-face in figure^ 

The four vertices have critical coordinates C^^ ={(0,0)}, C^2 = {(1,0), (0,1)}, 
C^3 = {(2,0), (1,2)}, Ct,4 = {(3,0), (0,1)}. The critical coordinates for the 
five edges areC^^^, ={(0,2), (2, 0)}, C^g^^ = {(2, 0)}, C^,^^ ={(0,2), (3,0)}, 
C-vxvA = {(0,2), (1, 1)}, C„3t,4 = {(3,0)}. There is just one 2— face with criti- 
cal coordinate (2, 2). The boundary operators of the simplicial chain complex 
associated to ^(3,2) o.f^ represented by the matrices 
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Figure 2: A non one-critical multifiltration 

and the matrices of the maps d are 
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3.1 One-critical case 

Following [3], a multifiltration such that there exists a unique critical coor- 
dinate for each a S B„(t>') and for all n is called one-critical. The module 
Cn is free if and only if the corresponding multifiltration is one-critical. In 
this case, 

Cn ^ fi'-^"' ^ -R"*" (3.3) 

and the boundary operator 9„ can be represented as a matrix with coef- 
ficients in R. It is also true that in the one-critical case the modules Z„ 
and Bn are naturally submodules of the free module Cn, thus it is possible 
to compute the quotient Hn = Zn/Bn in C„. Exploiting these properties, 
Carlsson et al. in [4J present an algorithm to compute a Grobner basis for 
multipersistence homology modules in the one-critical case. In general there 
is not a natural choice for an ambient free module, but by shifting Zn and 
Bn in Dn it is possible to compute the quotient x" Hn — Hn in Dn- Us- 
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ing this device, in section 5 we will generalize Carlsson's algorithm to all 
multifiltrations. 



Remark 2. In the one cviticdl case '. Cri 



Dn, is an isomorphism. 



Therefore, in this case, we have an identification of dn and dn up to isomor- 
phisms, as one can check from diagram 3.2, This gives the exact relationship 
between our approach and the one developed in 



4 One-critical case algorithm 

In this section we will recall the basic steps of the algorithm for the one- 
critical case [1]. For background information on Grobner basis theory we 
refer the reader to [TJ [9] . 

Let be a finitely generated free module with canonical basis {ei, ... cat}. 
Each element / G is a linear combination of monomials x'^ ej. The mod- 
ule generated by elements {fi, ... ft} in R'^ will be denoted hy < fi, ... ft >. 
We fix a monomial order |71|9] on the monomials of R^ . The leading mono- 
mial LM(/) and leading coefficient LC(/) of / € R^ are respectively the 
greatest monomial of / and its corresponding coefficient. If F is a finitely 
generated submodule of R^ , then LM.{F) is the submodule of R'^ generated 
by the leading monomials of the elements of F. 

Definition 4.1. A finite set of generators {/i, . . . , ft} for a module F C R^ 
is a Grobner basis of F if 

LM((/i, . . . , /t)) = (LM(/i), . . . , LM(/0) 

Using Buchberger algorithm, shown in listing [l} a Grobner basis of F C 
R^ can be computd from any finite set of generators. 

Definition 4.2. For each set of elements {fi,---,fm} £ , the syzygy 
module Syz{fi, . . . , fm) is the kernel of the map jp: R"" ^ R^ sending the 
canonical basis elements of R^ to {/i, . . . , fm}- 

A Grobner basis for the sygyzy module Syz{fi, . . . , fm) exists by Schreyer's 
theorem (see page 334 in [9]). We can compute such basis using the Wall 
algorithm described in [T^] or the Schreyer algorithm described in [1]. 

Buchberger and Wall (or Schreyer) algorithms share most of their logic. 
The Grobner bases for (/i, . . . , fm) and Syz{fi, . . . , fm) can indeed be com- 
puted concurrently as shown in listing [2] (using the Wall algorithm) . 

The algorithm described by Carlsson et al. in [4J uses these algorithms 
to compute Grobner bases for all Bn, Zn and Hn in the one-critical case. 
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The input of the algorithm are the matrices of the maps dn- The steps of 
the algorithm are the following: 

1. Compute the (reduced) Grobner basis of Bn using Buchberger algo- 
rithm on the columns of dn-^i- 

2. Compute the (reduced) Grobner basis of Zn as the syzygy module of 
the elements corresponding to the columns of dn- 

3. Compute the quotient Hn = Bn/Zn using the multivariate division 
algorithm. 

The Grobner bases of Bn and Zn+i can be computed concurrently using the 
algorithm we have described above. 

Remark 3. As stated in the previous section, the differential operators dn 
are matrices if and only if the multifiltration is one-critical. The previous 
algorithm therefore cannot be directly applied to all multifiltrations. 

In [Ij , Carlsson et al. observed that the entries of the matrices of dn are 
terms whose monomials only depend on the position in the matrix. They 
called the matrices with this property homogeneous, highlighting the fact 
that the image of dn is given by homogeneus elements of Cn-i- 

5 General Multipersistence Algorithm 

With our algorithm we compute, starting from general multifiltrations, re- 
duced Grobner bases of x" Bn, x" Zn and x" Hn — Hn for all n G {0 ... d}. 
To compute a Grobner basis for x" Bn we can use the Buchberger algorithm 
(as in step 1 of [4j ) on the columns of dn+i- 

To compute a Grobner basis for xf Zn we first compute the syzygy module 
of the elements corresponding to the columns of the matrix dn (as in step 2 
of [5j ) . This is however a submodule of R^^ . To obtain a Grobner basis of 
X Zn as desired, it is thus necessary to map the result in Dn- 
The multivariate dision algorithm (step 3 of [4j) can be used without any 
modifications to compute a Grobner basis of the quotient x"" Znjx" Bn- 

Remark 4. Observe that the matrices dn are homogeonous matrices, this 
implies that the computational complexity of our algorithm is polynomial- 
The algorithm in IjJ is in fact polynomial beacuse the matrices are homoge- 
neous and we perform the same operations, adding only multiplication by a 
matrix that is polynomial. 
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Example 5.1. We will now compute Grobner bases for the significant mod- 
ules arising from the multifiltration in figure^ In this example the only non 
trivial homologies are Hq and Hi . The modules Hq and Hi are isomorphic 
to Zq and Zi respectively. The information given by the couple {Zn, B^) 
is anyway more informative than the quotient H^ to understand the birth 
and death of generators in terms of multi-degree. As already stated, it is 
necessary to embed Zn and Bn in Dn to compute Grobner bases. 



The reduced Grobner bases for x^y'^Zo and x^y'^Bo are 
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The reduced Grobner bases for x y Zi and x y Bi are 
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Algorithm 1 Buchberger algorithm 
function Buchberger(/i, . . . , /„) 

Q ^ {/l, . • . SVECTOR(/,,/,) ^ Q \ Q < i < j < s} 

for all p G P do 

\i h = Reduce (p, G) / then 

Q ^gyj {h}, r^pu{ svector(5, h) ^o\geg} 

end if 
end for 

return Q 
end function 
function Svector(/i, /2) 

c ^ LC(/i)/LC(/2), Sij = lcm(LM(/i),LM(/,))/LM(/,) 

return S21 /i - CS12 /2 
end function 

function Reduce(/, {51, ... , gt}) 

while there exists a gi such that LM(^j) divides LM(/) do 
c ^ (LC(/) LM(/))/(LC(50 LM(50), f^f-cgi 

end while 

return / 
end function 
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Algorithm 2 Buchberger algorithm with syzygy computation 
function BuchbergerWithSyzygy(/i, . . . , f^) 

g^{{fi,ei),...,{fm,em)}, 5^0 
r ^ { SVECTOR((/i, €i), {fj,ej)) ^ I < i < J < s} 
for all p eV do 
P^V-{p} 

if {h, s) = Reduce(p, G) = (0, s) then 

S^Su{s} 
else 

g ^gu {h}, r^ru{ svector(5, h) ^o\geg} 

end if 
end for 

Q'^{g\{9,s)eg} 

return {g,S) 
end function 

function Svector((/i, si), (/2, .S2)) 

c ^ LC(/i)/LC(/2), Sij = lcm(LM(/i),LM(/,))/LM(/,) 

return (s2i fi - csi2/2,S2i si - csi2 S2) 
end function 

function Reduce((/, s),{{gi,si), {gt, st)}) 

while there exists a gi such that LM(g'j) divides LM(/) do 
(LC(/) LM(/))/(LC(50 LM(5,)) 
f^f-cgi, s-^s-csj 
end while 
return (/, s) 
end function 
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